Astronomy & Astrophysics manuscript no. CAbate 


©ESO 2013 


February 20, 2013 





Wind Roche-lobe overflow: 
Application to carbon-enhanced metal-poor stars 

C. Abated O. R. Pols', R. G. Izzard^, S. S. Mohamed^-^, and S. E. de Mink'^-^* 

' Department of Astrophysics/IMAPP, Radboud University Nijmegen, P.O. Box 9010, 6500 GL Nijmegen, The Netherlands 

e-mail: [CAbate; 0. Pols] Oastro. ru.nl 
~ Argelander Institut fiir Astronomie, Auf dem Hiigel 71, D-53121 Bonn, Germany 

e-mail: izzard@astro.uni-bonn.de 

South African Astronomical Observatory, Observatory Road, Observatory 7925, South Africa 

e-mail: mohamed@saao . ac . za 
* Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD, USA 

e-mail: deminkOstsci .edu 
' Dep. of Physics and Astronomy, Johns Hopkins University, 3400 N. Charles Street, Baltimore, MD 21218, USA 



Received accepted ... 

ABSTRACT 

Carbon-enhanced metal-poor stars (CEMP) are observed as a substantial fraction of the very metal-poor stars in the Galactic halo. 
Most CEMP stars are also enriched in .s-process elements and these are often found in binary systems. This suggests that the carbon 
enrichment is due to mass transfer in the past from an asymptotic giant branch (AGB) star on to a low-mass companion. Models 
of binary population synthesis are not able to reproduce the observed fraction of CEMP stars without invoking non-standard nu- 
cleosynthesis or a substantial change in the initial mass function. This is interpreted as evidence of missing physical ingredients in 
the models. Recent hydrodynamical simulations show that efficient wind mass transfer is possible in the case of the slow and dense 
winds typical of AGB stars through a mechanism called wind Roche-lobe overflow (WRLOF), which lies in between the canonical 
Bondi-Hoyle-Lyttleton (BHL) accretion and Roche-lobe overflow. WRLOF has an effect on the accretion efficiency of mass transfer 
and on the angular momentum lost by the binary system. The aim of this work is to understand the overall effect of WRLOF on the 
population of CEMP stars. To simulate populations of low-metallicity binaries we combined a synthetic nucleosynthesis model with 
a binary population synthesis code. In this code we implemented the WRLOF mechanism. We used the results of hydrodynamical 
simulations to model the effect of WRLOF on the accretion efficiency and we took the effect on the angular momentum loss into 
account by assuming a simple prescription. The combination of these two effects widens the range of systems that become CEMP 
stars towards longer initial orbital periods and lower mass secondary stars. As a consequence the number of CEMP stars predicted by 
our model increases by a factor 1.2-1.8 compared to earlier results that consider the BHL prescription. Moreover, higher enrichments 
of carbon are produced and the final orbital period distribution is shifted towards shorter periods. 
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1. Introduction 

The stellar population of the Galactic halo is characterised by 
low-mass stars with poor metal content. These stars are among 
the oldest that we observe, relics of the early stages of star for- 
mation in the MiUcy Way. The wide-field spectroscopic surveys 
HK ( |Beers et aL| ri992 ) and HES ( |Christheb et aLj[2gor| i, which 
are devoted to studying this population, reveal a surprisirigly 
high frequency of carbon-enhanced metal-poor stars ([C/FeVl> 
-I- 1.0, CEMP stars hereinafter) among the very metal-poor stars 
(VMP here indicating [Fe/H] < -2.0). The observed CEMP to 



VMP ratio is approximately 20% (for example: 25% Marsteller 
eTan [2007; 9 ± 2% jPrebel et all [20061, 21 + 2% [Lucatello 
et al.[ [2006 ), with the fraction of CEMP stars rapidly increas- 
ing for decreasing iron content and for increasing distance from 
the Galactic plane ( [CaroUo et aLl[2012| l. 



* Hubble fellow. 

' given two elements X and Y, their abundance ratio is [X/Y] = 
log|(|(A'x /A'y)-log|Q(A'xo/A'i'o)> where Nx,y refers to the number density 
of the elements X and Y and O denotes the abundance in the Sun. 



CEMP stars are classified into different groups according to 
the presence of the heavy elements barium and europium, which 
are produced by slow and rapid neutron-capture processes {s- 
process and r-process), respectively. The largest group of CEMP 
stars are the .s-process rich CEMP-,? stars, which display bar- 
ium enhancements of [Ba/Fe]> -1-0.5 and account for at least 



80% of all CEMP stars (Aoki et al. 2007 1. Among these stars 



some show enhancements of both r- and i-elements. A single 
case of a CEMP star highly enhanced only in r-process elements 
is also documented ( [Sneden et al. 2003 i. Finally, one group of 
stars does not exhibit pecuUar abundances of neutron-capture el- 
ements (e.g. [Aoki et al.] [2002| ). [Beers & Chris tHeb (2005) and 
'Masseron et al. (20101 provide detailed reviews of the CEMP 
subgroups. 

Several hypotheses have been put forward to explain the 
overabundances observed in CEMP stars: (a) the level of car- 
bon is primordial, or close to primordial, and was produced in 
the first generation of stars (e.g. Mackey et al.[ 2003 , and Cooke[ 



jet al. 201 1 1; (b) low-mass stars of extremely low metallicity 



might have undergone exotic mixing episodes that dredged up 
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internally produced carbon to the surface ( [Fujimoto et al.|[2000| l; 
(c) a binary scenario in which in the past carbon-rich material 
from a thermally-pulsing asymptotic giant branch (TPAGB) pri- 
mary star polluted the low-mass main-sequence secondary star 
Today we only observe the secondary. The primary has become 
an unseen white dwarf. The binary scenario is currently con- 
sidered the most likely formation mechanism for CEMP-s stars: 
in fact, the analysis performed by Lucatello et al.| ( [2005 b) on a 
sample of these objects demonstrates that the fraction of CEMP- 
s stars with detected radial-velocity variations is consistent with 
the hypothesis of all being members of binary systems. The same 
binary mass transfer scenario is invoked to explain the properties 
of Ba and CH stars (McClure & Woods worth 1990; ). 

A quantitative understanding of the origin of CEMP stars 
involves many branches of stellar astrophysics, some of which 
are still not well understood. The main uncertainties are related 
to (/) stellar evolution, particularly the nucleosynthesis during 
the AGE phase and internal mixing and diffusion processes in 
both stars of the binary system, (ii) the mass transfer process, 
(Hi) the binary fraction of low-mass stars in the Halo and the 
distribution of orbital parameters, (iv) the initial mass function 
(IMF) at low metallicity. 

Several studies have considered population models in an 
attempt to reproduce th e ob served CEMP fractio n of 9-25%. 

p.001) come to the 



Lucatello et al. (2005a) and Komiya et al. 



conclusion that an IMF biased toward intermediate-mass stars 
is required to reproduce the fraction of CEMP/VMP stars mea- 
sured in the Halo. However, large changes in the IMF are incon- 
sistent with the small fraction of nitrogen-enhanced metal-poor 
stars observed in the Halo (Pols et aT] |2012| l. With a different 
approach, Izzard et"ar| ( ,2009, 109 hereinafter) choose the solar 



neighbourhood IMF proposed by Kroupa et al. (1993 KTG93) 



in their population synthesis model and investigate the proper- 
ties of metal-poor stars. 109 try to reproduce the observed frac- 
tion of CEMP/VMF[^ stars by only varying uncertain physical 
parameters related to nucleosynthesis, mass transfer and mix- 
ing processes. For standard values of these physical parameters 
109 find a CEMP/VMP fraction lower than the observed one by 
almost a factor 10. With a set of parameters that reduces the 
minimum core mass required for third dredge up and allows ef- 
ficient third dredge up in stars of mass down to 0.8 Mq 109 find 
a CEMP/VMP fraction of approximately 9%, approaching the 
range of the observations. However, the observed distributions 
of carbon, nitrogen and heavy elements in CEMP stars are not 
well reproduced by the model. This suggests that the above men- 
tioned uncertainties need to be further investigated. In this paper 
we study the mass transfer process in some detail. 

Recently, different indications have emerged suggesting that 
the efficiency of wind mass transfer in binary systems has so 
far been underestimated, at least in the case of slow and dense 
winds characteristic of AGB stars. Wind mass transfer plays an 
important role in the binary scenario for the formation of CEMP 
stars, and, more generally, is involved in several problems of as- 
trophysical interest, e.g. the formation of CH and barium stars, 
the shaping of planetary nebulae, symbiotic stars, novae and the 
evolutionary path leading to the progenitors of Type la super- 
novae. Therefore in this paper we push forward the analysis of 
109 with a more accurate description of the wind mass transfer 
process and we investigate the effects of wind mass transfer on a 
population of binary stars, focusing in particular on CEMP stars. 

In Sect.l2|we briefly describe the context of wind mass trans- 
fer. In Sect.|3]we discuss the main parameters of our binary pop- 



ulation synthesis model and how wind mass transfer is imple- 
mented in our model. In Sect. |4]and[5] the results of our analysis 
are shown and discussed while Sect.lSlconcludes. 



2. Wind mass transfer 

In the binary scenario, the primary star produces carbon and s- 
elements during the late stages of its evolution, when it under- 
goes thermal pulses. Dredge-up episodes bring these elements 
to the surface of the star, where they may be lost in a wind or 
by Roche-lobe overflow, polluting the main-sequence compan- 
ion. The material accreted by the secondary star might in turn be 
diluted and burnt. When the donor is an AGB star, Roche-lobe 

T965) and will lead 



overflow is typically unstable (Paczynski 



the binary system to a common envelope phase, with negligible 
accretion to the secondary (see'Ricker & Taam 2008' for more 



details on accretion in a common envelope). Therefore the wind 
mass transfer scenario plays a crucial role in CEMP formation. 

2.1. Limitations of tlie standard scenario 



The canonical description by Bondi & Hoyle ( 1944 1 of the wind 
mass transfer mechanism is appropriate under the assumption 
that the wind velocity is much higher than the orbital velocity 
of the accreting star. This condition is not always fulfilled by 
AGB winds. In binary stars of periods around 10"* days the or- 
bital velocity is about 10 kms"', whereas outflows from AGB 
stars are observed with wind velocities in the range 5-30 km s~^ . 
The mechanism that drives the wind is not fully understood. It is 
thought that in the outer atmosphere of the AGB star, where the 
temperature is low enough, dust grains form and are accelerated 
by radiation pressure, dragging the surrounding gas along due to 
collisional momentum transfer ( |Hofne r 2009, Bladh & Hofner| 
2012). Therefore the ratio between the orbital velocity and the 
wind velocity, which in turn determines the applicability of the 
Bondi-Hoyle-Lyttleton (BHL) prescription, depends on the sep- 
aration and the mass of the two stars during the phase of mass 
transfer and on the chemical composition of the grains. 

2.2. The Wind Roche-iobe overflow mechanism 

In the past few years the possible existence of a relatively ef- 
ficient mode of wind mass accretion in a binary system has 
emerged both from observations and hydrodynamical simula- 
tions^ MiraAB ( [Karovska et al.j [2005] l and SS Leporis ( |Bhnd| 
jet al. 201 l| l are two examples of detached binary systems that 
are undergoing very efficient mass transfer (about the efficiency 
of mass transfer in Mka AB there is some debate, see e.g. 



ISokoloski & Bildsten 2010 1 although the AGB donor star is un- 
ambiguously filling only a fraction of its Roche lobe (approxi- 
mately 10% in the case of Mira, 80% for SS Lep). Moreover, 
recent hydrodynamical simulations suggest a new mode of mass 
transfer that lies somewhere in between Roche-lobe overflow 
and wind mass transfer. In this mode, called "wind Roche- 
lobe overflow" (|Mohamed & PodsiadlowsE 



2007 hereinafter 



In 109 stars with [Fe/H] < -2 are called EMP stars. 



WRLOF) or "gravitational focusing" ( |de Val-Borro et al.[[20(j9) , 
the wind of the primary star is focused towards the orbital plane 
and in particular towards the secondary star. This has two ma- 
jor effects on the evolution of a binary system. The first effect is 
that the accretion rate predicted in the WRLOF regime is signifi- 
cantly higher than the BHL predictions. The second effect is that 
most of the material that is not accreted is lost from the vicinity 
of the outer Lagrangian points L2 and L3, with consequences for 
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Fig. 1. Schematic picture of the WRLOF mechanism: R^, the 
wind acceleration radius, lies close to the Roche-lobe ra- 
dius. Inside the wind acceleration zone (shaded area) the wind is 
slow and can be efficiently accreted to the secondary through the 
first lagrangian point Li (sizes are not in scale). 



the angular momentum lost by the binary system that we further 
discuss in Sect. 13.3.41 

WRLOF can occur in systems where the wind is gravitation- 
ally confined to the Roche lobe of the primary star (the donor) 
and then falls into the potential well of the secondary through the 
inner Lagrangian point Li ( [Mohamed & Podsiadlowski] 2007| l. 
A good indicator of the conditions for the occurrence of WRLOF 
is the ratio Rd/RL.i, where: 



- ^L,i is the Roche-lobe radius of the primary star For a given 
mass ratio this quantity is proportional to the binary separa- 
tion [19^. 

- Rfi is the radius of the wind acceleration zone, namely the 
region where the wind is accelerated beyond the escape ve- 
locity. 

WRLOF can occur in systems where the wind acceleration ra- 
dius is larger than, or is a significant fraction of, the Roche-lobe 
radius of the wind-losing star. This condition is schematically 
represented in Fig.[T| where the wind acceleration zone is repre- 
sented as a shaded area around the star of radius 

In AGB stars we assume the acceleration of the wind is 
driven by radiation pressure on dust grains and therefore /Jj co- 
incides with the dust formation radius. In this case R^ scales with 
the stellar radius, and it also depends on the effective tempera- 
ture of the AGB star, Tf-g, and on the chemical composition of 



the dust. Hofner (2007i suggests the following approximate re- 
lation: 

^ \ cond / 

Tmnd is the condensation temperature of the dust, a parameter 
that depends on the chemical composition of the compound, for 
example: Tcond = 1500 K for carbon-rich dust (C/0> 1) and 
T^cond = 1000 K for oxygen-rich dust (C/0< 1). During most of 
the AGB phase the dust formation radius is linearly proportional 
to the stellar radius and for carbon-rich dust Rd/R* ~ 3. 

[Mohamed (2010) uses a smoothed-particle hydrodynamics 
(SPH) code to simulate wind mass transfer in Mira-like circu- 
lar binary systems with a IM© primary in the AGB phase and 



a 0.6 Mq secondary with separations between 5 and 60 AU. A 
similar approach is used by |de Val-Borro et al. ( 2009) but with 
Ml = 1.2 Mq, M2 = 0.6 Mo and a varying up to 70 AU. The 
evolution of each binary is followed in detail at least until the 
wind material has expanded to twice the size of the primary's 
Roche-lobe radius, and typically for an amount of time equal to 
many orbital periods. This amount of time is considered long 
enough to study the development of wind anisotropics, spiral 
shocks, equatorial flows, and to follow the fate of the lost ma- 
terial, distinguishing the fraction of mass that is accreted from 
the fraction that is lost or remains bound to the system but is not 
accreted. 

The wind is modelled by inserting particles in shells at the 
surface of the star, where the time interval between successive 
injections is regulated to produce a continuous outflow of mate- 
rial and to minimise discontinuities between the shells. Different 
types of wind are modelled: the one that best reproduces the 
properties of an AGB star like Mira considers a wind outflow 
with initial speed of 5.5 kms ' at the stellar surface, and in- 
cludes radial pulsations implemented using a piston approxima- 
tion (e.g. Freytag & Hofner 2008) and radiation pressure on dust 
grains. In every model the evolution in time of the positions and 
fluid properties is followed, predicting among other parameters 
the dust-formation radius, the ratios of the density, velocity, wind 
anisotropy at the equator relative to the poles, and the accretion 
efficiency. 

Several parameters are likely to influence the mass transfer 
mechanism in Mira-type binaries, first of all the mass-loss rate, 
which is determined by the luminosity, temperature and pulsa- 
tion period of the wind-losing star, although it is chosen to be 
approximately 10"''MQyr"' in the simulations. The model for 
the formation of the dust grains and the C/O ratio affects the 
velocity of the winds and therefore the intensity of WRLOF fea- 
tures. Moreover, stellar rotation and magnetic fields are likely to 
play an important role in shaping the outflow. 

In Table [T| we list orbital periods, Ra/Rla ratios and accre- 
tion efficiencies ySacc, defined as the ratio between the mass ac- 
creted by the secondary star and the mass lost by donor star per 
unit of time. These values are the result of the simulations of 
Mo hamed] ( |2010| l with the model of the AGB wind described 
earlier in this Section in which the evolution of every binary sys- 
tem is followed for at least one orbital period. Tests performed 
by Mohamed] f 20 1 0) with different assumptions about the wind 
driving mechanism give an estimate of the uncertainty in /3.^cc 
which is within 50% of each value (e.g. for P - 1.34 x 10^ days 
/3acc =0.10 ±0.05). 

The data in Table [T] are shown as plus signs in Fig. |2] to- 
gether with a proposed model (solid line) which is further dis- 
cussed in Sect. |3.3[ High Rd/Rh.i values correspond to close 
systems whereas low values correspond to wide systems be- 
cause Rd/RL,i <^ Rt/a- From wide to close systems (i.e. from 
left to right in Fig. |2]i the accretion efficiency ySacc initially in- 
creases, then reaches a maximum and finally decreases again 
with Rd/RL,\ - This behaviour is explained as follows: for a given 
Ml and M2, in wide systems the wind acceleration zone is 
smaller than the Roche lobe, therefore the wind is fast and only 
a small fraction of it is accreted, as also predicted by the BHL 
prescription. Moving to smaller separations (i.e. greater /?d//?L,i) 
the dust formation radius becomes a progressively more signif- 
icant fraction of the Roche-lobe radius, the wind is increasingly 
confined in the primary's Roche lobe and therefore the accre- 
tion efficiency grows. With smaller a the accretion efficiency de- 
creases because in close systems a large fraction of the wind es- 
capes through L2 and Lj and is not accreted to the secondary star. 
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Table 1. Initial period P, ratio between the dust formation ra- 
dius and the Roche-lobe radius, Ra/Rla, and accretion effi- 
ciency ySacc for binary systems with a 1 .0 Mq AGB primary and 
a 0.6 Mq main-sequence companion (from [Moham ed 2010i). 
Because a (and does not depend on a), high values 

of the Rd/RL,\ ratio correspond to close systems and low values 
correspond to wide systems. 



P/days 




^acc 


1.34 X 10' 


0.40 


0.10 


5.98 X lO'* 


0.67 


0.24 


2.58 X lO'* 


1.18 


0.45 


9.12 X 10^ 


2.4 


0.35 


3.21 X 10' 


2.8 


0.10 



We refer to |Mohamed & Podsiadlowski| ( |2007l l and to |Mohamed| 
( |2010| l for further details about the WRLOF mechanism. For 
the purposes of our work we apply the results of the above- 
mentioned hydrodynamical simulations to our binary population 
synthesis code, as described later in Sect. 3.3 and we evaluate 



how this mode of mass transfer affects the predictions for CEMP 
stars in a population of binary systems. 

3. Models 

In order to perform this study we make use of the binary popula 
tion synthesis code described by |Izzard et al.| ( |2004||20()6||2009 



Our model follows binary evolution according to the rapid bi- 
nary stellar evolution prescription of Hurley et al. ( 2002[l and 
the algorithms for AGB evolution and nucleosynthesis of |Izzard| 
|et al . (2004, 2006). For a complete discussion of our model, its 
characteristics, parameters and uncertainties we refer to 109. In 
this section we briefly summarise the most important charac- 
teristics of our model (Sect. 3.1 1, the updates that we introduce 



compared to the work by 109 (Sect. 3.2 1 and the way we imple- 
ment WRLOF (Sect. 133). 



3.1. Population synthesis 

Our population synthesis simulations are based on a grid of A^^ 
binary evolution models distributed uniformly in In Mi - In M2 - 
In a parameter space, where Mi and M2 are the initial masses 
of the primary and of the secondary star respectively, a is the 
initial separation of the system and we take - 128. The ini- 
tial metallicity of our model sets is Z = 10 or equivalently 
[Fe/H]= -2.3. We consider circular orbits, therefore the eccen- 
tricity, e, is always zero. 

We count the stars of a particular type, for example CEMP 
stars, according to the sum 



^type 



Ml.min M2,mi„ 



6M1 6M2 6a 5t , (2) 



"min 'imn 



where, as in the paper by 109: 



- the size of a cell in the parameter space is 5M\ ■ 6M2 ■ 5a; the 
timestep is 6t. 

- 5 is the star formation rate which is assumed to be constant; 

- <5type is equal to 1 when the star is of the required type and 
is zero otherwise. Stars are selected from our model popula- 
tion according to their age, surface gravity and surface abun- 
dances, as follows: VMP stars are older than 10 Gyr and 
with surface gravity logio(g/cms"^) < 4.0 ; CEMP stars are 



VMP stars characterised by a surface abundance of carbon 
[C/Fe]> 1.0; CEMP-s stars are CEMP stars also enriched 
in barium, [Ba/Fe]> 0.5. 

- Ml and M2 vary respectively in the ranges [0.7, 8.0] M© and 
[0.1,0.9] Mq. Initially M2 < Mi by definition. 

- a varies between 3 and 10^ Rq. To be able to compare our 
work with the results of 109 we assume that all stars are 
formed in binary systems with this range of separations. In 
reality some VMP stars are single or in wider orbits and we 
can take these stars into account by reducing the binary frac- 
tion. 

- t varies in the range [10, 13.7] Gyr, the approximate ages of 
the Halo and the Universe. 

- *!* is the inital distribution of Mi, M2, and a. We assume that 
*!* is separable. 



= >I'(Mi,M2,fl) = ip{M,) mDxia) 



(3) 



where the primary mass distribution i/r(Mi ) is the intial mass 
function by KTG93, the secondary mass distribution <piM2) 
is flat in <7 = M2/M1 (any mass ratio < ^ < 1 is equally 
likely), the separation distribution ;^(a) is flat in In a (i.e. 
Xia) cc I /a). 

Our model sets assume efficient thermohaline mixing: the ac- 
creted material mixes instantaneously with the stellar envelope. 
This approximation is reasonable in many cases, as suggested 
by the calculations of Stanchffe et al. ( 2007|l, eve n though more 
recent studies (e.g. [Stancliffe & Glebbeek 2008[ ) show that the 
situation is more complicated and other processes such as gravi- 
tational settling in some cases reduce the effect of thermohaline 
mixing. 



3.2. Parameter ciioices and updates 

The nucleosynthesis algorithm which follows the evolution of 
the star through the first, second and third dredge ups, modi- 
fying the surface abundances as required, is mostly based on 
the work by Karakas et al. ( 2002[ l and Karakas & Lattanzio 



( |2007| l. A prescription for hot-bottom burning is also i ncluded 



Abundances are scaled according to the solar values by I Anders 



& Grevesse ( 1989|l. Win d mass-loss rates are parameterised ac 



cording to the |Reimers| ( |1975 1 formula multiplied by a factor 



T] = 0.5 on the first giant branch andlVassiliadis & Wood 
on the AGB. The '^C pocket efficiency f 13, as defined by 



T993]l 



zzard 



|et al.| ( [200 6^, is set to 1 by default. The common envelope evo 



lution is treated according to the prescription of [Hurley et al. 



(2002 ) with acE - 1, and we do not include accretion during the 
common envelope phase. Our default model set A uses the same 
input physics as the model set A described in Sect. 2.1.5 and 4.1 
of the paper by 109. We hst here the updates we introduce. 

- Wind velocity. To calculate the mass accretion rate, 109 
assume the wind velocity is a fraction of the escape velocity 
from the mass-losing star. However, the escape velocity 
decreases as the radius of the star and its mass-loss rate 
increase, whereas the wind velocities observed in TPAGB 
stars show an opposite trend. Therefore to calculate the mass 
accretion rate when the primary is a TPAGB star we take the 
wind velocity defined in Eq. (3) of the article by Vassiliadis] 
& Wood ( 1993| ), which was derived by interpolating the 
observations. We choose a minimum value of Vw = 5 km s~\ 
which is consistent with the observations (e.g. |Knapp &| 
|Morris|19831 ). As [Vassiliadis & Wood, (1 993 j we introduce a 
maximum value of Vw - 15 km s"\ 
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Angular momentum loss. The algorithm used by 109 to cal- 
culate the change to the orbital angular momentum is based 
on Eq. (21) of the article by Hurley et al. (2002). This equa- 
tion only applies if the mass lost by the wind is much greater 
than the mass accreted by the secondary star. In most of our 
model sets we calculate the change to the orbital angular mo- 
mentum assuming a spherically symmetric wind, 



[[Mr 



•M, 



)m2 + (M2w 



■M, 



'orb 



(Ml + M2)2 



(4) 

where Qoib - 2jt/P is the orbital angular frequency and 
M,w and M,a are, respectively, the wind mass-loss rate and 
accretion rate of star /. When the accretion efficiency is low, 
as in the BHL prescription, Eq. (21) of Hurley et al. (2002 1 
approaches Eq. (j4]). Therefore we do not expect the CEMP 
to VMP ratio to change significantly after this modification 
of the model. In fact, the fraction of CEMP to VMP stars 
that we compute with the two equations is the same: 2.30% 
and 2.31% (using respectively Eq. (21) of Hurley et al.| 
2002 , and Eq. (HV, see Table |2|. The model sets A' - Cq and 
(J - Jq in Table]/] include Eq. ^ and we will refer to these 
as iiw-model sets. 




Rjj/Rrl 

Fig. 2. Models of the wind accretion efficiency ySacc for different 
values of q calculated with Eq. (]9]). The solid line is for q = 
0.6, the mass ratio used by Mohamed (2010). The dot-dashed 
and dotted lines show the resulting efficiency for q = Q.l and 
q = I, respectively. Plus signs are the results of hydrodynamical 
simulations by Mohamed (2010[) Usted in Table]!] 



3.3. A proposed model for Wind Roche-lobe Overflow 

In the WRLOF mode of mass transfer the wind material is 
ejected in a highly aspherical geometry. This has two major ef- 
fects, on the accretion efficiency of the mass transfer and on the 
angular momentum lost by the binary systems. We discuss our 
prescription to consider these two effects in Sect. 3.3. 1 and Sect. 
3.3.4 respectively. 



3.3.1. The WRLOF accretion efficiency 

To use the results of the hydrodynamical simulations in our 
model, we fitted the data of Table [T] with a parabolic function. 



Cl X + C2 X + C2 



(5) 



where x = Ra/Rla, cj = -0.284 (+0.018), C2 = 0.918 (±0.057) 
and C3 = -0.234 (±0.034). The errors result from the non-linear 
least-squares algorithm used for the fit and do not convey the 
uncertainty associated to the simulations. Since in none of the 
models calculated by Moha med] ( |2010 ) the accretion efficiency 
is higher than 50% we impose a maximum value ySacc.max - 0.5 
on ySacc- The result of the best fit is shown in Fig.]2]as a solid line. 

Since Eq. (]5]l can assume zero or negative values in the range 
of possible values of Rd/Rui, we always choose the maximum 
accretion rate between WRLOF and BHL. This implies that for 
large separations, i.e. Rd/Rh,! smaller than about 0.4, we calcu- 
late the accretion rate with the formulation of the BHL prescrip- 
tion given by [Boffin & Jorissen| ( |1988| l, 



BHL 



2vr^ 



GM7 



1 + 



Voib 



(6) 



where a - 1.5 is a constant, Vorb is the relative orbital veloc- 
ity, G is the gravitational constant and the other quantities have 
the meaning already described. This is consistent with Mohamed 
(2010 1, in which at large separations the same accretion rates are 
obtained as BHL. 



0.5 



0.4 



0.3 



0.2 



BHL 
WRLOF 




10' 

Initial Period/days 

Fig. 3. Efficiency of wind mass transfer ySacc as a function of the 
initial period for a binary system with Mj = 1 .0 Mg and M2 = 
0.6 Mq calculated with BHL (solid Hne) and WRLOF (dotted 
line). 



3.3.2. An example of accretion efficiency in tine WRLOF 
regime 

To illustrate the effects of WRLOF on the wind mass trans- 
fer process we compute the accretion efficiency in binary sys- 
tems with different initial periods and fixed initial masses. Mi = 
1.0 Mg and M2 = 0.6 Mg as in the simulations of |Mohamed| 
( 2010[ l. In Fig. 3]we plot the accretion efficiency calculated with 
the BHL (solid line) and the WRLOF (dotted line) models as a 
function of the initial period. The maximum accretion efficiency 
in the case of WRLOF is almost 50% and this high mass transfer 
efficiency is possible over a wide range of periods. A 1 Mg star 
loses approximately 0.2 Mg during the AGB phase, therefore its 
companion can reach 0.7 Mg, close to the minimum mass for a 
metal-poor star to evolve off the main sequence in a Hubble time. 
The maximum accretion efficiency reached in the BHL model is 
roughly 15% for periods around 1800 days, which means the 
companion accretes at most 0.03 Mg and typically much less. 
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Table 2. Physical parameters in our binary population models and the derived CEMP/VMP ratio. Model set A is the same as in 
the article by 109. All models use the wind velocity prescription of VassiUadis & Wood ( 1993 i except model sets A and A' where 
Vw I'esc- The errors convey only Poisson statistics and could be reduced by increasing the resolution of our grid of models. A 
finer grid is not necessary because the observed CEMP/VMP ratio varies between 9% and 25%, a much wider range than the errors 
in our models. 



Model set 



Ang. mom. loss 



Third dredge up CEMPA^MP (%) 



A as Hurl ey et al. ( 2002) BHL 

A' ^~ssw - bq. lEJ BHL 



standard 
standard 



2.30: 

2.31 : 



0.04 
0.04 



B 
C 

D 

E 

G 
J 

I 

K 

K„ 



ssw - Eq 
ssw - Eq. I? I 
ssw - Eq. I?! 
72-Eq. 1^1 
72 -Eq. iTOl 
72-Eq. iTTJl 
ssw - Eq. (?i 
ssw - Eq. I?! 
ssw - Eq. (5 1 
72-Eq. 1^1 
72-Eq. iTUi 
72-Eq. iTOi 



BHL 
WRLOF 
WRLOF g-dep 
BHL 
WRLOF 
WRLOF g-dep 
BHL 
WRLOF 
WRLOF ^-dep 
BHL 
WRLOF 
WRLOF ^-dep 



standard 
standard 
standard 
standard 
standard 
standard 
enhanced 
enhanced 
enhanced 
enhanced 
enhanced 
enhanced 



2.22 4 
2.90 4 
2.63 4 
3.124 
4.06 4 
3.85 4 
9.00 4 

11.15 : 
11.91 : 
10.27 : 
12.96: 
13.66: 



0.04 
0.04 
0.04 
0.04 
0.03 
0.04 
OJT" 

:0.13 
:0.14 
:0.13 
:0.14 
:0.14 



3.3.3. Mass-ratio dependence of the WRLOF efficiency 

Eq. (j5]l implicitly assumes that the accretion efficiency ySacc is 
independent of the mass ratio q - M2/M1. No conclusion can 
be drawn about any dependence on q based on the results of 
Mohamed ( 2010 1 which are calculated for a fixed mass ratio q - 
0.6. Therefore, additional assumptions are required to introduce 
a relation between jSacc and q. In Eq. (j6]l the dependence on q for 
a fixed Mi can be expressed as. 



2vr 



GM, 



GM, 



(7) 



Hence the accretion efficiency /3bhl scales with a factor q^ and 
there is an additional dependence related to the term T in the 
square brackets. If term T dominates, i.e. for small separations 
when Vorb > Vw, Eq. (j7]| becomes: 



a 



GMi 



2^fl^ {l+qfl^ \av: 



1/2 



(8) 



Vice versa when T is not important, e.g. for large separations, 
then ySfiML scales simply as the mass ratio squared, q^. For sim- 
plicity we choose to consider only the latter dependence in the 
model of WRLOF, keeping in mind that large separations corre- 
spond to small orbital velocities and therefore to the case when 
Voib ^ v^, the assumption under which the BHL prescription is 
appropriate. 

Assuming this dependence on q of the mass-accretion effi- 
ciency in the WRLOF we introduce the following expression: 



min 



25 



q \c\ X + C2 x + c-i 



j 5 Pace, ma: 



(9) 



where x, ci, C2, C3, ySacc.max are as in Eq. dSll and the factor 25/9 
arises because Eq. (j9]) should coincide with Eq. (j5]l for q - 0.6. 
In Fig.|2]we show our models for different values of q. 

We point out that the term T is not negligible in general. 
Therefore we expect the real dependence offlacc on the mass ra- 
tio to be weaker than q^. We consider Eqs. (jsj) and ^ to bracket 
the dependence on q of the accretion efficiency. 



The calculations of Mohamed ( 2010[ l are for AGB primaries 
with oxygen-rich winds and dust, i.e. C/O < I. Therefore the 
value of the parameters Ci, C2 and C3 of the fit could be differ- 
ent in the case of carbon-rich winds (C/O > 1), which is the 
relevant situation for CEMP stars. However, at the moment we 
cannot take this into account because there are no hydrodynam- 
ical simulations of binary systems with a carbon-rich primary 
star. 

3.3.4. Angular momentum loss in the WRLOF regime 



The hydrodynamical simulations performed by Mohamed 
( |2010 ) show that the density of the matter lost in the wind is 
strongly enhanced towards the orbital plane. This enhancement 
is clear also in wider systems where the accretion efficiency is 
consistent with the BHL predictions. Because the geometry of 
the wind is not spherically symmetric, a prescription alternative 
to Eq. Q is required to calculate the amount of angular momen- 
tum carried away by the wind. 

Let us assume that the material lost from the binary system 
carries away a multiple y of the average specific orbital angular 
momentum. 



J - yx 



Jo 



Ml + M2 



(Mi + M2) . 



(10) 



If the matter is lost from the vicinity of L2 and L3 one expects 
efficient loss of angular momentum with 7 > L However, the 
exact amount of angular momentum lost in a WRLOF situa- 
tion is not well constrained by the hydrodynamical models of 
Mohamedj (j20T0|. [Jahanara et al.| (]2005| have studied angular 
momentum loss associated with wind mass transfer, also using 
hydrodynamical simulations but with different assumptions re- 
garding the wind mechanism. In the case of a 'radiative wind' 
- which corresponds best to the case of an AGB wind - with 
a velocity much lower than the orbital velocity, they find that 
matter leaving the binary system has an average specific angular 
momentum of approximately 0.6 x a^Qorb- This corresponds to 
a value y x 2.4-3.2 for mass ratios in the range 1/3-3. In their 
study of the formation mechanism of barium stars Izzar d et al.| 
(2010 ) use Eq. (10 1 with y = 2 to reproduce the periods and ec- 
centricities of the barium stars. We use the same prescription in 
some of our models, i.e. in model sets D - Eq, I, K, Kq in Table|2] 
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and we will refer to these as 72 -model sets. In these models we 
apply efficient angular momentum loss with y = 2 to all binary 
systems because the geometry of the wind is modified also in 
systems that are not in the regime of enhanced accretion. 

We note that angular momentum loss from a spherically 
symmetric wind from the primary star, Eq. (jij) with M2W = 0, 
can be written in the form of Eq. (lOi with y - M2/M1. 



Compared to this ssw situation, our y2 models with y = 2 corre- 
spond to a mode in which the binary system loses more orbital 
angular momentum and therefore shrinks more. 



4. Results 

We start our analysis with model set A which has the same input 
physics as the default model A in the article by 109 and predicts 
the same CEMP to VMP ratio of 2.30%. We successively modify 
model set A with the changes mentioned in Sect. 3.2 and 3.3 



and we study the effect of each of these updates on the overall 
population of CEMP stars. In Table |2] our different model sets 
are listed with the corresponding CEMP to VMP ratio. 

4.1. Model set B: wind velocity 



The wind-velocity relation of Vassiliadis & Wood ( 1993 1 im- 
plies that v^v increases with time, reaches a maximum of 15 
kms ' in the regime of superwind and then remains constant 
until the end of the AGB phase. On the other hand, in model set 
A the wind velocity is directly proportional to the escape veloc- 
ity, I'esc, which decreases in time. Even though the dependence 
of the wind accretion efficiency on Vw is strong (see Eq.[6]), in the 
regime of strong mass loss the actual difference between the two 
calculated values of is small, within a factor 2. Therefore we 
find that in model set B the CEMP to VMP ratio is only slightly 
reduced to 2.22%. 

4.2. Model set D: efficient angular momentum loss 



The hydrodynamical simulations of |Mohamed (2010 1 show that 
the geometry of the wind ejected by the primary star is modified 
by the gravitational field of the companion in a wide range of 
systems, including those that are not in the regime of enhanced 
accretion efficiency. It is therefore instructive to study the effect 
of efficient angular momentum loss on a population of CEMP 
stars independently. In model set D the variations in the orbital 
angular momentum are computed using y = 2 in Eq. ( 10 1, while 
the accretion rate is computed with the BHL prescription Eq. (j6|. 
Model set D predicts a CEMP to VMP ratio of 3.12% and com- 
pared to model sets A and A' produces notable differences in the 
distribution of initial and final period of the CEMP population. 
The initial period range of systems that lead to a CEMP star is 
broadened and shifted towards wider separations (see Fig. [4] top 
panel). Systems shrink more due to wind mass loss and there- 
fore wide binaries become close enough to efficiently transfer 
material. On the other hand, initially close binaries tend to un- 
dergo a common envelope phase, in some cases merge, and do 
not become CEMP stars. Analogously the final period distribu- 
tion is shifted towards relatively close systems and the distribu- 
tion peaks at a period more than one order of magnitude shorter 
than in set A, as shown in the bottom panel of Fig. [4] 

A direct comparison of our population with a measured pe- 
riod distribution is not possible: the number of CEMP stars with 
known period is low, despite the fact that many of them are 
known to be binaries, because long periods are difficult to mea- 



< 
z 

< 




10' 
Period/days 



Fig. 4. Initial (top panel) and final (bottom panel) period distri- 
butions of the CEMP population. In model set D (dotted line) an- 
gular momentum losses are computed with Eq. ( 10 1 and y - 2. 



Solid line is the standard model set A. Plus signs are the peri- 
ods of observed stars. The bins are equally spaced and the width 
in log[0(P/days) of each bin is 0.1 (top panel) and 0.2 (bottom 
panel). In this and in the plots that follow the y-axis indicates 
the expected number of CEMP stars in each bin and the plots 
are normalised such that the area under the graph is the same for 
each model and model A peaks at 1 . 



sure. Since the average time for which surveys have been ongo- 
ing is about ten years, this is approximately the maximum pe- 
riod that we can presently measure. In the observational SAGA 
database only six CEMP stars have a measured period (see Table 
4 in the article by 109), and they are displayed in Fig. |4]as plus 
signs (two systems have almost the same period, 1.26 x lO"* and 
1.29 X 10^ respectively, and therefore are indistinguishable in 
the plot). The two shortest-period CEMP binaries are not repro- 
duced by any of our models. The characteristics of the distribu- 
tions of other parameters, e.g. Mi, M2, [C/Fe], are not signifi- 
cantly different from model set A. 

4.3. Model sets with WRLOF accretion efficiency 

In this section we present the results that we obtain after adopt- 
ing our model of accretion efficiency in the WRLOF mode of 
mass transfer. If we consider a spherically symmetric wind from 
the donor star, the CEMP to VMP ratio is 2.63% when the 
WRLOF accretion efficiency is proportional to q~ (model set 
Cq), and 2.90% in the case of ^-independent WRLOF (model set 
C). The corresponding 7'2-model sets, dubbed Eq and E, predict 
a CEMP to VMP ratio of 3.85% and 4.06% respectively. 
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Fig. 5. a) Distribution of initial secondaiy mass in binaries that lead to the formation of a CEMP star, b) Distribution of initial period 
of binaries that form CEMP stars. Angular momentum variations are calculated with the spherically simmetiic wind assumption 
(Eq. |4]i. The 0.1 dex bins are equally spaced in M2/M0 (left panel) or logio(P/days) (right panel), c) and d) As in (a) and (b), we 
„ compute angular momentum loss. Model sets Cq and Eq (dot-dashed lines) are the result of ^-dependent 

!q. (|9]l. Model sets C and E (dotted lines) are the result of ^-independent WRLOF. The solid lines show our 



use y = 2 in Eq. ( 10 



WRLOF defined by 
default model set A. 



4.3.1 . Model sets C and Cqi CEMP initial parameter space 

The distribution of initial Mi is essentially the same as the cor- 
responding distribution shown in Fig. 4a of 109: it tells us that 
the majority of CEMP stars are formed by wind accretion from 
a primary of 1.2 - 1.5 Mg. Below this limit the primary does not 
undergo third dredge up, and therefore no carbon is brought to 
the surface and transferred to the companion. The highest pro- 
duction of carbon during the TPAGB phase occurs in the mass- 
range 2 - 2.5 M0, but these stars are much less numerous, ac- 
cording to the IMF of KTG93. 

Most notable is the distribution of initial secondary mass 
shown in Fig. |5^: each of the three model sets peak at 0.8 M© 
but C and Cq show a considerable fraction of CEMP stars that 
come from systems with very low M2 (dotted and dot-dashed 
line respectively). This low-mass tail is absent if we consider 
the BHL mode of accretion (model set A, solid line), because 
with the WRLOF prescription the secondary can accrete mass 
more efficiently. For example a 1 .5 M© primary expels around 
0.8 Mq of carbon-rich material during its TPAGB phase: hence, 
accreting 50% of the lost matter even a star of about 0.4 Mg in- 
creases in mass sufliciently to evolve to logm(g/cms"^) < 4.0 
when f > 10 Gyr and be selected as CEMP star If we observe 



today such a CEMP star, a large amount of mass we see is from 
the AGB primary and not the original low-mass companion. 

The difference in shape between model sets C and Cq re- 
sults from the mass-ratio dependence of model set Cq. In this 
set, for a given primary mass a low-mass secondary star accretes 
a smaller amount of material than a more massive secondary star. 
Therefore the low-mass tail of model set Cq does not extend to 
less than 0.42 Mg, whereas in set C it extends down to 0.15 Mg. 

Examining the distributions of initial periods in Fig.|5]3, the 
effect of WRLOF is to increase the number of wide systems that 
lead to CEMP star formation. WRLOF accretion is efficient in 
a wider range of separation than BHL, as illustrated in Fig. [3] 
Compared to set Cq, in model set C more CEMP stars come from 
systems with initial period between 10,000 and 50,000 days. 
This difference is due to the secondary stars coming from the 
low-mass tail. To become a CEMP star, a low-mass star needs to 
accrete material from a relatively massive primary with a large 
dust formation radius. Therefore such a binary system under- 
goes the conditions for WRLOF only when the initial separation 
is relatively wide. 
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Fig. 6. Distribution of [C/Fe] in the CEMP population. The his- 
togram shows the observed distribution in our data sample, with 
Poisson errors. The solid, dotted and dot-dashed lines are model 
sets A, Cq and Eq respectively. Model set A is our default, Cq 
and Eq are the models for ^-dependent WRLOF. Model set Cq 



assumes a spherically symmetric wind. Model set Eq 
lated with 7 = 2 in Eq. (fTO]). 



4.3.2. Model sets E and Eq: CEMP initial parameter space 

The distributions of M2 resulting from model sets E and Eq (Fig. 
[5]:) are qualitatively similar to those obtained from sets C and 
Cq (Fig.|5^). However, the low-mass tails of the distributions of 
M2 do not extend to as low masses as in the previous sets. A less 
massive secondary star needs to accrete more mass in order to 
become a CEMP star, hence it is more likely to form in a sys- 
tem with a relatively high mass primary. Such systems are more 
likely to fill their Roche lobe and undergo a common-envelope 
phase in our y2-model sets because of stronger angular momen- 
tum loss. Therefore CEMP stars forming from very low-mass 
secondaries are relatively rare in the 'y2-model sets. 

The strong angular momentum loss causes the formation of 
a large fraction of CEMP stars from very wide systems, as can 
be seen in the initial period distribution (Fig. 5]l). This effect was 
already appreciable in the case of BHL wind mass transfer (Fig. 
[4]) and with the WRLOF mode it is further increased. 

The comparison of the initial parameter distributions allows 
us to compare the effects of the WRLOF and BHL mechanisms, 
although initial masses and periods are not observable. In the 
next sections we study the distributions of carbon and nitro- 
gen that are directly comparable to the distributions observed 
in CEMP stars. 



4.4. Distributions of carbon and nitrogen 

To compare our models with observations we use a database of 
VMP stars based on the SAGA observational database compiled 



by Suda et al. 



(2008 



2011 1, combined with data of metal-poor 



stars from |Frebel et ai.| ( |20b6l ) and from |Lucatello et aL] ( |2006| l. 

We select stars corresponding to our VMP criteria: (/) the iron 
abundance is in the range [Fe/H] - -2.3 + 0.5 dex and (//) the 
star is a giant or a sub-giant, with logjQ(g/cm s"^) < 4.0 (in the 
database of |Frebel et al. 2006 stars are catalogued as giants or 
sub-giants according to their B - V colour). This selection leaves 
us with a database of 413 VMP stars with measured carbon abun- 
dance, of which 100 are CEMP stars. Our database collects all 
data of VMP and CEMP stars so far available in the literature 
and is not complete because of the different observational prop- 



erties and selection effects of the original samples: e.g. data from 
Lucatello et al. (2006i contains mostly CEMP stars. Therefore 
the CEMP fraction in our database (approximately 24%) is not 
representative. 

Our mass-ratio dependent models Cq and Eq predict abun- 
dance distributions of carbon and nitrogen very similar to our 
^'-independent models C and E respectively. Therefore in this 
section we focus our discussion on the model sets Cq and Eq. 

None of our model sets is able to accurately reproduce the 
distribution of carbon in the observed CEMP population, as 
shown in Fig.l6] Our model sets Cq and Eq show qualitatively the 
same trend as the data. However, the peak around [C/Fe]^ 2.5 
corresponds to a carbon abundance 3-4 times higher than the 
observations. This might indicate that in our models the accre- 
tion efficiency is too high. The fate of the transfered material 
is to be mixed by the thermohaline process in the envelope of 
the accretor, composed mainly of hydrogen. However, a star of 
0.5 Mo must accrete 0.3 - 0.4 Mq from the donor to be selected 
as a CEMP star, and this material is only weakly diluted in the 
envelope. Hence secondary stars with low initial mass become 

we estimated the er- 



is calcu- strongly carbon enriched stars. In Sect. 2.2 



ror on the accretion efficiency calculated in the hydrodynami- 
cal simulations to be within 50%. If we reduce the efficiency 
of WRLOF by a factor 2 in our model set Cq we remove the 
stars with initial mass below 0.6 M© from the CEMP population 
and the CEMP/VMP ratio is consequently reduced from 2.63% 
to 2.53%. This also results in a shift of the distribution towards 
lower values of [C/Fe] but the peak remains around 2.1 - 2.3, 
i.e. higher than the observations. A further decrease in the ac- 
cretion efficiency in our WRLOF model would presumably re- 
sult in a better agreement with the observed carbon distribution. 
However, we should keep in mind that the data in our observed 
sample may not be representative of the real carbon distribution 
because our sample is inhomogeneous and incomplete. 

In Fig. |7] we show the distribution of [N/Fe] versus [C/Fe] 
for y2-model set Eq. Most of our CEMP stars have small nitro- 
gen enhancements and we predict a small number of CNEMP 
stars, i.e. CEMP stars enhanced in nitrogen ([N/Fe] > 1.0 and 
[N/C] > 0.5). Most observed CEMP stars are in the region be- 
tween these two groups, where we predict a dearth of systems, 
and no CNEMP stars are found in our data sample. Our pre- 
dictions are similar to those of 109 and this is not unexpected, 
because our updates modify the efficiency of the mass transfer, 
not the chemical composition of the transferred material. To re- 
produce the nitrogen enhancement some extra mixing processes 
are required to burn part of the carbon in the CN cycle and mod- 
ify the amount of nitrogen brought to the surface of the AGE star 
2003] l. 



(see e.g. NoUett et al. 



4.5. Model sets witli efficient tliird dredge-up at low mass 

In our model populations the fraction of CEMP stars is consid- 
erably lower than that observed. Our failure to reproduce the ob- 
served distributions of carbon and nitrogen might simply be a 
consequence of this. To increase the number of CEMP stars pre- 
dicted by their standard model set A, 109 choose a set of param- 
eters that reduces the minimum core mass and envelope mass 
required for third dredge up and increases the amount of mate- 
rial dredged up. With such a set of parameters all TPAGB stars 
down to initial masses of 0.8 M© undergo efficient third dredge 
up and 109 find a CEMP/VMP ratio of 9.43%, about four times 
higher than in their standard model A and at the low end of the 
range of the observations. 
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3.0 



Fig. 7. The distribution of [N/Fe] vs [C/Fe] in the CEMP popu- 
lation predicted by y2 -model set Eq with ^-dependent WRLOF. 
Darker grey indicates greater density of stars. Crosses are data 
from our selection of SAGA database. Average error bars of the 
observations are shown in the top-left corner CEMP stars appear 
to the right of the vertical dotted line. CNEMP stars enhanced 
both in carbon and nitrogen appear above the diagonal dotted 
line. 



oT 1.5 



O 




Fig. 8. As Fig. |6] for models with enhanced third dredge up. 
Model G (solid line) is calculated with the BHL prescription, 
Jq (dotted line) and Kq (dot-dashed line) are ^-dependent model 
sets. 



The need for efficient dredge up at low core masses and 
low metallicity has been pointed out by many authors, e.g. 
Groenewegen & de Jong ( 1993); Marigo et al. ( 1996) and |Izzard| 
& Tout,(2004f . Further justification for these assumptions is pro- 



vided by recent detailed models that find dredge up for stel- 



lar mass down to 0.9 M© at Z = 10 (see e.g. Stancliffe & 



[Glebbeek, 2008 Karakasl|201 0|l. Indirect indications of efficient 
dredge up at low mass and metallicity arise from the observed 
abundance of s-elements in metal-poor stars ( Bonacic Marinovic 
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Fig. 9. Period distribution in our CEMP population. The 0.2 dex 
bins are equally spaced in log [q(P/ days). Solid, dotted and dot- 
dashed line are model sets A, Eq and Kq respectively. Plus signs 
are the observations from SAGA database. A is the default set. 
Eq and Kq are 72-model sets with <7-dependent WRLOF. Model 
set Kq assumes efficient dredge up in stars down to 0.8 Mq. 



|et al.| |2007| l and from observations of low-mass white dwarfs in 
the globular cluster M4 (see [Kalirai et al.||2009] l. 

In our model set G we take ffie same set of parameters as in 
model set G of 109 in addition to the equation for the wind veloc- 
ity proposed by VassiUadis & Wood ( 1993] ) and the assumption 
of a spherically symmetric wind and we obtain a CEMP/VMP 
ratio of 9.00%. The small difference with respect to the result 
of 109 is due to the different equation for the wind velocity we 
use to calculate the accretion rate. When we include the WRLOF 
prescription in our model G we obtain the same relative increase 
as in the models without enhanced dredge up. The WRLOF ac- 
cretion efficiency increases the CEMP star fraction by a factor of 
about 1 .2- 1 .3, which further increases to approximately 1 .4- 1 .5 
when we also consider the contribution of efficient angular mo- 
mentum loss. With the ^-dependent WRLOF prescription the 
CEMP/VMP ratio increases to 11.91% in the ssw-v(\odt\ set Jq 
and to 13.66% in the y2-model set Kq. If we assume a binary 
fraction of 100%, these results are within the range of the ob- 
servations (9 - 25%). However, our models predict most of the 
CEMP stars to have a carbon-enhancement [C/Fe] > 2.5, much 
larger than the observations (Fig.|8]l. Because this choice of pa- 
rameters does not directly affect the amount of nitrogen trans- 
ported to the stellar surface the discrepancy in nitrogen abun- 
dances remains an issue. 

4.6. Orbital periods 

In Fig. [9] we compare the final period distributions from our y2- 
model sets Eq and Kq to our default model A and to the data 
from our selection of the SAGA database. Model set Eq predicts 
most CEMP stars at orbital periods around 2000 days, whereas 
default model set A predicts mostly wider systems and peaks 
around 40,000 days. Model set Kq, which allows efficient third 
dredge up in low-mass stars, predicts CEMP stars in a range of 
orbital periods wider than model set Eq and is able to reproduce 
the shortest-period CEMP stars of our sample. A similar result 
is obtained with model sets G and Jq. 

Even though many CEMP stars are known to be binaries, 
only a few have known orbital periods because long periods are 
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difficult to measure. In our observed sample 6 out of 100 CEMP 
stars show a period below 4400 days (the longest-period mea- 
sured in a CEMP star so far), i.e. a fraction of 0.06 if we assume 
that all the other observed CEMP stars are in binary systems 
with a longer period. In reality this is a lower limit because only 
a fraction of the observed systems in our heterogeneous sample 
have been investigated to determine the period. Model sets A, Eq 
and Kq predict a fraction of CEMP stars with period shorter than 
4400 days of 0.2 , 0.48 and 0.36 respectively (values in the same 
range are found with all the other model sets). This exceeds the 
observational lower limit, and therefore we are not able to rule 
out any of our models. 

5. Discussion 



This paper presents the first study in which the WRLOF mode of 
mass transfer is taken into account in binary population synthe- 
sis: the purpose of this work is to investigate how the implemen- 
tation of a more realistic description of the wind mass transfer 
process affects a population of CEMP stars. 

WRLOF has two effects on the evolution of binary systems: 
it increases both the accretion efficiency of the mass transfer and 
the angular momentum taken away by the material lost in the 
wind. The enhanced accretion efficiency produces an increase in 
the CEMP to VMP number ratio by a factor 1.2 - 1.3. It has 
strong effects on the initial secondary mass distribution, which 
gains low-mass secondary stars, and on the carbon distribution, 
which peaks at large [C/Fe]. Our prescription of efficient an- 
gular momentum loss increases the CEMP to VMP ratio by 
a factor 1.2 - 1.5 and has strong effects on the initial and fi- 
nal period distributions of CEMP stars. Initially wide systems 
shrink enough to interact and transfer mass to the secondary star 
Moreover, a higher number of systems are predicted at relatively 
short final periods. The combined contribution of these two ef- 
fects widens the initial range of systems that become CEMP stars 
towards longer periods and lower secondary masses. As a result 
the CEMP to VMP ratio predicted by our models increases by a 
factor 1.4 - 1.8 and we find a distribution of carbon abundances 
more similar to observations than in previous studies. 

The reliability of our model of WRLOF is limited by the 
small number of hydrodynamical simulations of wind accretion 
in binary systems. In particular, the efficiency of angular mo- 
mentum loss by the stellar wind is poorly constrained by cur- 
rently available models. We have assumed two different, very 
simple models: a spherically symmetric wind that does not in- 
teract with the binary system (equivalent to y = ^ < 1 in eq. 10 1, 
and a strongly modified wind that carries away twice the average 
specific angular momentum of the orbit (y = 2). These assump- 
tions may respresent very wide systems and relatively close sys- 
tems, respectively, because the specific angular momentum car- 
ried away by the wind is likely to decrease wiffi increasing ra- 
tio of wind velocity over orbital velocity (Jahanara et al.||2005|. 
The dependence of y on orbital period and mass ratio cannot 
be quantified further at present, and our simple assumptions can 
only bracket the real situation. To improve upon our model of the 
enhanced accretion rate and of the efficiency of angular momen- 
tum loss more hydrodynamical simulations are needed, and for 
a wider range of initial masses and separations. Considering the 
uncertainties related to the accretion rate, we tested how strongly 
the CEMP to VMP fraction depends on the accretion efficiency 
in the WRLOF regime. The dependence is found to be weak: if 
the accretion efficiency is reduced by a factor 2, the CEMP/VMP 
ratio decreases from 2.90% to 2.77% in model set C and from 
2.63% to 2.53% in model set Cn. 



Another uncertainty is related to the condensation tempera- 
ture used to determine the dust formation radius of AGB stars. 
A variation in Tcond affects the range of separations over which 
binary systems undergo WRLOF. We tested this by varying 
Tcond by 10% around the assumed value of 1500 K. A lower 
condensation temperature (Tcond = 1350K) does not affect the 
CEMP/VMP ratio while the initial and final period distributions 
include more wide systems. If T^ond is 1650 K the CEMP/VMP 
in model C is decreased from 2.90% to 2.50% because the num- 
ber of CEMP stars with initial period longer than 30,000 days is 
almost zero and initially close systems, with period shorter than 
1000 days, enter into a common envelope which prevents mass 
accretion. Our results are also affected by unceitainties related 
to the AGB phase, such as the amount of carbon and nitrogen 
dredged up to the surface and the wind acceleration mechanism, 
to the mixing processes that occur in the primary and the sec- 
ondary star and to the IMF at low metallicity. 

We have assumed that all stars are formed in binary systems 
with semi-major axes 3 < a/R© < 10^. The observed fraction 
of binary systems with this range of orbits in ffie solar neigh- 
bourhood is smaller, however. |Duquennoy & Mayor ( 1991 1 find 
that in the solar neighbourhood about 60% of solar-Uke stars 
have a binary companion. The study by |Kouwenhoven et"aL] 
( 2QQ7\ of ffie young stellar association Scorpius OB2 shows a 
binary fraction of at least 70% (3cr confidence) and probably 
close to unity among stars of spectral type A and B in the range 
5Rq < fl < 5 X IO^Rq. With the assumed flat distribution in 
log JO fl, the range of separations in our grid corresponds to about 
75% of the range analysed in Scorpius OB2. This implies that 
our calculated CEMP/VMP ratios should be reduced by a fac- 
tor < 0.75, although we note that the fraction and distiibution 
of binary systems among metal-poor halo stars is not well con- 
strained and may be different than in the solar neighbourhood. 

WRLOF can influence a variety of astrophysical phenomena 



outside the context of this paper (Mohamed & Podsiadlowski 



2007). For instance, WRLOF could likely play a role in the 
wind-accretion scenario which is generally accepted to explain 
the formation of the barium and CH stars observed at relatively 
higher metallicity. Also, in the case of formation of planetary 
nebulae in binary systems WRLOF widens the range of sepa- 
rations in which the presence of a companion star distorts the 
shape of the ejected material. Finally, in a symbiotic binary the 
white dwarf could accrete enough mass from an AGB compan- 
ion star to trigger nova outbursts or the explosion of a Type la 
supernova. It would be interesting to study the effects on these 
phenomena of our implementation of WRLOF. 



6. Summary and conclusions 

We summarise here our analysis. 

We have updated our population synthesis code by improv- 
ing algorithms to calculate the angular momentum varia- 
tions and wind velocity. These changes result in a modest 
change in the overall population of carbon-enhanced metal- 
poor (CEMP) stars. 

//) We propose a model for the wind Roche-lobe overflow 
(WRLOF) mode of mass transfer. In this model the accretion 
efficiency of mass transfer is based on the results of hydrody- 



namical simulations of Mohamed (20101 and only depends 



on the radius of the donor star and on the separation of the 
binary system. We also introduce a dependency on the mass 
ratio of the two stars that scales as q^. We consider the effect 
of WRLOF on the orbital angular momentum lost by the bi- 
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nary systems with a simple prescription in which the wind 
carries away twice the average specific angular momentum 
of the orbit. 

Hi) The WRLOF accretion efficiency causes an increase by a 
factor 1.2-1.3 of the number ratio of CEMP stai's to VMP 
stars. The prescription for efficient angular momentum loss 
increases the CEMP/VMP ratio by a factor 1.2-1.5. When 
both effects are taken into account we predict a CEMP/VMP 
ratio up to 4% if we assume a binary fraction of unity. 

iv) Our model with efficient third dredge-up in stars of mass 
down to 0.8 Mq predicts a CEMP/VMP ratio of 9%, con- 
firming the results of jlzzard et a l. (2009). When we com- 
bine efficient dredge up with the WRLOF model we find 
a similar relative increase in the CEMP/VMP ratio as pre- 
dicted in the models with standard dredge up and we predict 
a CEMP/VMP ratio in the range 1 1 - 14% (again for a binary 
fraction of unity). 

In conclusion, our model of WRLOF widens the range of 
binary systems which form CEMP stars and as a result in- 
creases the predicted CEMP fraction and modifies the distribu- 
tion of carbon abundance and period of the CEMP population. 
Assuming in addition efficient dredge up in low-mass stars we 
obtain CEMP/VMP ratios in the range of the observations. Other 
physical processes need to be considered to reproduce the distri- 
bution of carbon and nitrogen observed in CEMP stars. 
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